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A traditional view on solitons in optical fibers as robust particle-like structures suited for informa- 
tion transmission has been significantly altered and broadened over the past decade, when solitons 
have been found to play the major role in generation of octave broad supercontinuum spectra in 
photonic-crystal and other types of optical fibers. This remarkable spectral broadening is achieved 
through complex processes of dispersive radiation being scattered from, emitted and transformed 
by solitons. Thus solitons have emerged as the major players in nonlinear frequency conversion in 
optical fibers. Unexpected analogies of these processes have been found with dynamics of ultracold 
atoms and ocean waves. This colloquium focuses on recent understanding and new insights into 
physics of soliton-radiation interaction and supercontinuum generation. 
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I. INTRODUCTION 

Interplay of dispersion and nonlinear self-action in 
wave dynamics has been at the focus of attention across 
many branches of physics since the middle of the past 
century after the seminal Fermi-Pasta-Ulam work on heat 
dissipation in solids. This work has been quickly followed 
by the discovery of many nonlinear wave equations inte- 
grable with the inverse scattering technique (1ST. Soli- 
tons is a particular class of localized and remarkably ro- 
bust solutions found with the 1ST technique. Soliton 
studies have quickly become a subject of its own and have 
soon developed far beyond the initial subset of integrable 
models. Localised nonlinear waves in the non-integrable 
models are often called solitary waves to distinguish them 
from 1ST solitons. However, the use of the original term 
'soliton' has now spread widely into non-integrable cases. 



FIG. 1 (color online) Photograph of supercontin- 
uum generated by a 800nm femtosecond pulse in a 
photonic crystal fiber and projected onto a screen. 
(http://www.bath.ac.uk/physics/groups/cppm/). 



Solitons have been experimentally observed and stud- 
ied theoretically in fluid dynamics, mechanical systems, 
condensed matter and notably in nonlinear optics. His- 
toric and scientific accounts of these developments can be 
found, e.g., in the excellent book by Alwyn Scott (|ScottL 

Op t ical solitons in fibers (jHasegawa and TappertL 
1 19731 : iMollenauer et all I198Q[ ). have been most ex- 
tensively researched as potential information carriers 
(JAgrawai I2QQ7I : IMollenauer and Gordonl . l2QQ6h . With 



this view in mind solitons can be treated as particle-like 
objects and their dynamics can be conveniently reduced 
to the Newton-hke equations for the soliton degrees_of 
freedom, such as, e.g. position and phase 



20071: iGorshkov and Ostrovskyi [l98l[ iKaup and Newelj 
19781 : IMollenauer and Gordonl . l2006l ). Ability of per- 



turbed solitons to emit dispersive radiation is a prop- 
erty vividly expressing their wave nature. Radiation 
emissi on by solitons has been well known in 20th cen- 
turv ([Akhmediev and Karlssonl Il995 : iKarpmanl . 1 19931 ; 
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iKivshar and Malomedl . Il989l : IWai et a/.Ul990l ). however, 
it was considered at that time as something that to the 
large degree undermines usefulness of the soliton con- 
cept. This vision has changed dramatically around and 
after 2000, when the first experimental observation of 
an octave wide spectral broadening in photonic crysta l 
fibers was reported by Ranka et al. (|Ranka et all |2000[ ) 
and subsequently reproduced in dozens of labs. This 
effect has become known as the fiber supercontinuum , 
see, e.g.. (iDudlev e t al!. '2006' 'Pudlev a nd Tavlorl 120091 : 
iKnight and Skrvalpin . 2007 : Russell , 2006[ ) for recent re- 
views, and Fig. [H The above reviews discuss impressive 
apphcations of supercontinuum for frequency comb gen- 
eration, in metrology, spectroscopy and imaging. 

Already after first experiments on generation of super- 
continuum in fibers it has become obvious that solitons 



are the major players in this process. Spectrahy they 
dominate in the infrared, where the group velocity dis- 
persion (GVD) is typically anomalous. The visible part 
of supercontinua most often spans through the range of 
normal GVD and is associated with dispersive radiation. 
Crucially, the dispersive radiation and solitons overlap in 
the time domain and hence interact by means of the Kerr 
and Raman nonlinear it ies. 

What was overlooked in the research on fiber soli- 
tons in the pre-supercontinuum era is the fact that 
interaction of solitons with dispersive waves and the 
associated frequency conversion processes can be ef- 
ficient and practically important. Not only experi- 
ments, but also theory of the soliton-radiation inter- 
action was not developed much beyond the spectrally 
narrow results following from the integrable Nonlin- 
ear Schrodinger (NLS) equation, where the solitons 
are largely inse nsitive to the inter a ction with other 
wave s (Kivshar and Maloniedl, [l989l : iKuznetsov et all 
I1995I : IZakharov and ShabatI , Il972[ ). Our primary aim 
here is to explain the frequency conversion effects result- 
ing from the soliton-radiation interaction and leading to 
supercontinuum generation in optical fibers, and to dis- 
cuss few other fascinating soliton related effects, which 
have stem from the supercontinuum research. 



II. MODELING OF SUPERCONTINUUM 

Talking about supercontinuum generation and soliton- 
radiation interaction we consider a fiber with sufficiently 
high nonlinearity and the zero GVD point close to the 
pump wavelength. These conditions are found in fibers 
with silica cores of few (~ 1 — 5) microns in diameter, 
which can be, e.g., either photonic crystal or tapered 
fibers pumped with a variety of sources. The later include 
femtosecond pulses of mode-locked Ti:sapphire lasers 
with wavelength around 800nm (Ranka et al.^ 2000) and 
nano-second pulses from microchip lasers close to 1/im 
(jStone and Knightl , |2008[ ). An important aspect of the 
dispersion profile allowing to achieve wide supercontinua 
is the normal GVD extending towards shorter ('bluer') 
wavelengths, see Fig. [2fa). Spectral broad ening with the 
opposite GVD slope has al so been studied (JEfimov et all 
l2QQ4HHarbold et a/J . [2002l ). 

The widely accepted model well reproducing experi- 
mental measurements of supercontinuum is 



a) 



Wavelength (|a.m) 
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d,A = ik{idt)A + (1) 

POO 

ij{l - 0)\A\^A + ijOA / dt'R{t')\A{t - t')\^. 
Jo 

Here z is the propagation length and t is the time mea- 
sured in the reference frame moving with the group ve- 
locity at the pump frequency. 7 ~ 0.01/W/m is the 
nonlinear parameter. = 0.18 measures the strength 
of the Raman nonlinearity with respect to the Kerr one 
and R(t) is the Rama n response function (Agr awali l2007t 
iGorbach et alluOOw . A is the modal amplitude, which 
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FIG. 2 (color online) (a) Group index, that is, speed of light 
divided by the group velocity, Ug = cdsk, (dashed line) and 
GVD P2 = dsk (full line) typical for photonic-crystal fibres 
used in supercontinuum experiments. The zero GVD point 
is at ^ 790nm (380THz) and the normal GVD range is for 
wavelengths < 790nm. (b) XFROG spectrogram showing si- 
multaneous frequency and time domain pictures of the super- 
continuum from Fig. [3] at z = 1.5m. Pump wavelength is 
850nm (350 THz). 



carries all the spectrum we are interested in. k{idt) is 
the dispersion operator, usually represented as a poly- 
nomial of the order N in idt. Dispersion profile of a 
linear fiber is recovered by taking A = e^K^)^-'^^^ ^ where 
^(^) — Xln=2 ~\Pn^^ 1 which is a polynomial fit (not a 
spectrally local Taylor expansion) of the fiber dispersion. 
Positive ^'s correspond to the absolute frequencies larger 
than the pump frequency. GVD coefficient ^2 is calcu- 
lated as /32 = 9|/c. /32 = determines location of the zero 
GVD frequencies. 

Eq. ^ includes three ingredients: 

• dispersion (1st term in the righthand side) 

• nonlinear phase modulation due to instantaneous 
Kerr nonlinearity (2nd term) 

• Raman scattering (3rd term) 

Taken separately these effects are well known. Notably, 
the last two lead to generation of new spectral content 



FIG. 3 (color online) Numerical simulation of supercontin- 
uum generation in a photonic-crystal fibre pumped with 200- 
fs pulses at 850 nm and having 6-kW peak power. The fiber 
dispersion as in Fig. [2ja). One can see 3 stages in the super- 
continuum expansion. Symmetric spectral broadening due to 
SPM happens over first 5cm. Antisymmetric spectral broad- 
ening due to soliton fission accompanied by the emission of 
the resonant Cherenkov radiation and soliton-radiation inter- 
action both develop between ^ 5cm and ~ 15cm. After that 
the dynamics of the short wavelength edge is determined by 
the radiation trapping. 



acting on their own. However, conversion into discrete 
Raman side-bands does not manifest itself in the typi- 
cal supercontinuum experiments. Also self- (SPM) and 
cross- (XPM) phase modulations (as particular cases of 
the generic nonlinear phase modulation) generate spec- 
tra, which are relatively narrow (Agrawal, 2007), and 
can not explain fully developed supercontinua. It is only 
when dispersion comes into play and all three of the 
above effects act together in symphony, then the qual- 
itatively new phenomenon of supercontinuum generation 
occurs. A typical supercontinuum experiment with the 
femtosecond 800nm pump leads to spectra covering range 
from 400nm to 2000nm afte r propagation in a ^ Im . 
long photonic crystal fiber (.Gu et all l20Q3l : iRanka et all 
I2OQOI : IWadsworth et all 120021). Figs, [^b) and [3] show 
numerical modelling of this process using Eq. ([!]). 

Results obtained using Eq. ([1]) have been di- 
rectly compared with many experimental measurements 
of femtosecond pulse propagation in photonic crystal 
fibers in both time and frequency do mains and ex- 
cellent agr e emen t has been established (lEfimoy et al. , 
200^, '2006', '2005'; 'Go rbach et ^ . l2006l : iGu et all 12003 : 
Skryabin et al , 2003 ). Probably the most convenient 



2007bll3: iHill et all l2QQ9l: iKudlinski and Mussoti . l2008l : 



Stone and Knightl . boOSi : iTraversl \2m¥ . 



way of representing the data for such comparisons is us- 
ing the cross-correlation frequ ency resolved optical gating 
(XF ROG) spectrograms (Efi mov et all l2004l : \Gn et all 
[2OO3I : iHori et all l2004l ). The XFROG spectrogram is 
the fourier transform of a product of the signal field 
E with the reference pulse Er delayed by time r: 
Ixfrog{t, 5) = I J^^ Er{t - T)E{t)e-'''dt\\ Er is usu- 
ally a pump pulse and the product of the two fields is 
generated by the sum frequency process in a x^^^ crystal. 
As with any model, Eq. ([1]) has limits to its ap- 
plicability. In particular, it is not applicable to de- 
scribe sharp field variations happening over few fem- 
toseconds and less. However, so far the role of such 
ultrashort features in typical fiber supercontinuum ex- 
periments has not been revealed and the octave spec- 
tral broadening happens through nonlinear interactions 
of dispersive waves and solitons which both are well de- 
scribed by Eq. ([1]). In particular, formation of the 
coupled soliton-radiation states with the continuously 
blue shifting radiation component and the red shift- 
ing s oliton plays a key r ole in supercontinuum expan- 
sion ([Cumberland et aU l2008l : iGorbach and SkrvabinI , 



The relative smoothness of a temporal field profile as- 
sociated with a typical supercontinuum is o ne of the 
reaso ns why self-steepening (wave breaking) ([Agrawal 
I2OO7I ) does not make a notable quantitative and qualita- 
tive impact if included or excluded from Eq. ([!]). Fur- 
thermore, going beyond the amplitude model is likely 
to be ne eded for correct descriptio n of self-steepening, 
see, e.g., ("Amirana shvili et al|j2008l ). Therefore we have 
opted for not taking it into account. In our experience, 
the most noticeable disagreements of Eq. ([1]) with the 
experimental measurements are due to omission of the 
frequency dependent losses. Other neglected effects are 
possible excitation of higher order or orthogonally polar- 
ized modes, dispersion of the Kerr nonlinearity, third har- 
monic generation and noise. Spontaneous Raman noise 
and input noise are unavoidably present in experiments. 
The spectra averaged over the statistical ensemble may, 
in some cases, be better suited for a comparison with 
measurements, where the fine spectral features of a sin- 
gle pulse excitation can be obscured due to resolution and 
long (ove r many pulses ) integration time of a spectrum 
analyser ([Corwin et all l2003l : iDudlev and Coenl . l2002h . 
Small to moderate pulse to pulse fluctuations in the gen- 
erated spectra should not prevent us from unraveling and 
understanding of the deterministic nonlinear processes 
underlying the supercontinuum generation and the soli- 
ton role in it. 



III. SOLITON SELF-FREQUENCY SHIFT AND SOLITON 
DISPERSION 

During the first stage of the supercontinuum devel- 
opment SPM induces spectral broade ning and associ- 
ated chirping of a femtosecond pulse ([Agrawaj l2007l ), 
which are compensated by the anomalous GVD (9|/c < 
0). Competition of these processes initiates formation 
of multiple solitons from an initial high power pump 
pulse (soliton fission), see Fig. [3l Amongst the host 
of the soliton related effects, there are two, which are 
particularly important for us. These are the soliton 
self- frequency shift induced by the Raman scattering 
([Mitschke and Mollenauerl , Il986l ) and emission of disper- 
sive radiation. The latter comes from the overlap of the 
soliton spectrum with t he range of frequencies wfi e re the 
fiber GVD is normal ([Akhmediev and KarlssonI , 1 19951 : 
'Karpmam ^ 19931 : IWai et all Il990[ ). 

Spectra of femto- or pico-second solitons are suffi- 
ciently narrow, so that the Raman gain profile of silica 
peaking at 13THz can be approximated by a straight 
line rising across the soliton spectrum from the nega- 
tive values (damping) at the short- wavelength end of the 
soliton spectrum to the positive values (gain) at the long- 
wavelength end (AgrawaJ, l2007l : iGorbach and SkrvabinI , 
I2OO8I : iLuan et all l2006h . The net result on the soliton 
is that its spectral center of mass is shifted towards red- 



der frequencies (Mitschk e and Mollenauerl , 1 19861 ) , which 
is referred as the sohton self- frequency shift. The group 
index Ug (ratio of the vacuum speed of light c to the 
group velocity, Ug = cdsk) increases with the wavelength 
providing that the GVD is anomalous. This leads to con- 
tinuous negative acceleration (deceleration) of solitons by 
the Raman scattering, which is very important for super- 
continuum, but often forgotten in soliton studies. 

An approximate soliton solution of Eq. ([1]) mov- 
ing with a constant acceleration can be derived under 
the assumption that the fiber group index varies lin- 
early in frequency (i.e., dsk ^ 6 and hence k ~ 6'^). 
For the zero initial frequency, this solution is given 
bv (iGamon and Belangerl 1 19901 : iGorbach and Skrvabinl . 
l2007alR): 



As = il){t — ts) exp [—iSst - 
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The soliton delay caused by the Raman effect is ts = 
g^z^ j2. g' > is the soliton wavenumber shift propor- 
tional to its intensity. The soliton trajectory in the (t, z\ 
plane is a parabola given by t = t^ (see Fig. [H^a)). Jg 
is the soliton frequency, which decreases linearly with z 
(see Fig. IHb)). 

Soliton-radiation interaction discussed below is sensi- 
tive to the phase matching conditions. Therefore it is 
important to derive spectral representation for an accel- 
erating soliton. The soliton spectrum is calculated as 



~As{^) 



/CO 
diA.e*** = e*'^(^)+"^[''-''^ltA((5 - (5,), (6) 
-CX) 



where '0((^) = J^ ip{t)e^^^dt is a real function. 
Wavenumbers of the soliton spectral components are 

ks{S) = dM^)^[S-Ss]ts] = q^dsk\s=6A^-Ss]^k{Ss). 

(7) 
ks{S) is linear in J, expressing the fact that the soliton 
is immune to GVD and hence d'^ks = 0. Also it can be 
seen that ks (S) is actually a tangent to the dispersion of 
Hnear waves k{d) = |/32^^ taken dit 6 = Ss and shifted 
away from it by the offset q. If GVD is anomalous for all 
the frequencies {/32 < 0) the soliton spectrum does not 
touch the spectrum of linear dispersive waves, see Fig. 
He). 

?/;, as in Eq. (J3j), is of course an approximate solution, 
and the Raman effect forces the soliton to shake off some 
radiation, see Fig. [2fb), Fig. IH^a). This radiation can 
be approximated by Airy function s (JAkhmediev et all 
1 19961 : iGorbach and Skrvabinl . |2QQ8[ ) . Airy waves emitted 
in a typical supercontinuum setting are very weak, broad 
band and have spectrum almost exclusively belonging to 




FIG. 4 Soliton evolution in {z,t) (a) and {z,6) (b) planes in 
the presence of the Raman effect. k{S) = ^/32^^. Note that 
the dispersive waves in (a) are not reflected by the soliton. 
(c) shows phase matching diagrams between the soliton and 
dispersive waves (see Eq. ()1Q|) ). Straight lines are the soliton 
wavenumbers ks for two different frequencies Ss (see Eq. [7|. 
Point 1 indicates the initial soliton frequency 6s and point 2 
is Ss at some distance down the fiber. Parabola shows the 
dispersion of linear waves k{6). All units are dimensionless. 



the anomalous GVD range. The latter is the main rea- 
son why they practically do not interact with solitons 
and their impact on the supercontinuum spectrum is sec- 
ondary in importance. 



IV. RADIATION EMISSION BY SOLITONS 

Solitons dominate the long-wavelength edge of 
the supercontinuum (see Figs, lb, 2), while the 
so-ca lled resonant (o r 'Cherenkov') radiation 
( Akh mediev and KarlssonL 



1995 : iBiancalana et 



i2QQ4 .Karpmanl . ll993Tiskrvabin et all I2QQ3I : IWai et al 



1990D comes to mind as one of the reaso ns for the spec- 
truni created in the normal GVD ra nge (Cri stiani et all 
120041 : iGu et ali 120031 : iHerrmann et aL ^20021). Despite 
the importance of the resonant radiation identified 
i n the first efforts t o m odel fi ber supercontinuum 
(jHerrmann et all 120021 : iHusakou and Herrmannl . l200l[ ). 
later it has become clear that it is only when the 
initially emitted radiation has a chance to interact 
with the solitons over long propagation distance, the 
expanding supercontinua ob served in experimen t s can 
be reproduced in mod ell ing (iGentv et all [2004ai 12005 ; 
Gorbach and Skrvabinl . l2007bl la: ISkrvabin and Yulin 



2005[ ). The Raman effect is the key factor ensuring such 



interaction. 
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FIG. 5 The same as Fig. [2 but with /^s > 0. {k{6) = 
1/^2^^ + ^PsS^)- Intersections of the straight hnes with k{S), 
see (c), correspond to the Cherenkov resonances. 
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FIG. 6 The same as Fig. H but with /^s < 0. {k{6) = 
^p2S^ + |/33(5^). Intersections of the straight hnes with k(6), 
see (c), correspond to the Cherenkov resonances. 



To model Cherenkov radiation it is often sufficient to 
include third order dispersion only, so that a single zero 
of GVD is taken into account 



id,A 



Mt + T^P^dt ^ = 



2V^ t 31 



(8) 



Linear dispersive solution of Eq. (J8j) is A = e*^^ ^^^ with 
K^) = \P2^'^ + \M^' The zero GVD point is located 
at ^0 = — /^2//^3- The soliton existence condition is /32 < 
and hence if, e.g., j3^ > 0, then the GVD is normal 
(9|/c > 0) towards higher frequencies 5 > Sq^ see Figs. 
[Sjc). The pulse spectrum entering this range is not able 
to propagate in the soliton regime, which leads to the 
radiation emission. 

Small amplitude Cherenkov radiation Fch with 
frequency detuned far f rom th e soliton obeys 
(lAkhmediev and Karlsson , 1995 : iBiancalana et at. . 




iKarpmanl . Il993l : JSkrvabin et a/.Ll2QQ3l : IWai et all . 



id. - \m - ^/^3s,^ 



FcK ^ dlAs 



(9) 



Unlike the above mentioned Airy waves, the resonant 
radiation is phase matching dependent and therefore 
it is narrow band. Using Fourier expansion Fch — 
J^ ech{^)^^^^^^^~^^^ d5 ^ taking Eq. ([6|), and calculating 
z derivatives of the phases involved, we then equate the 
wave numbers of the left and right hand sides of Eq. (|9|) 
to find the wavenumber matching condition 



m = k,{8) 



(10) 



The roots of Eq. (p!Q|) are the resonance frequencies 
^Ch' /^3 > (/^a < 0) gives the blue (red) shifted dis- 
persive waves, see Figs. [51 [6l respectively. Accounting 
for the higher order dispersions is e ffortless and may 
lead to several roots of Eq. (ITQI) (iFalk et al\ . 



m- 



iFrosz et all l2QQ5l : iGentv et all l2QQ4b[ ). Energy of the 
emitted wave (or waves) is drawn from the entire soliton 
spectrum resulting in the adiabatic decay of t he soliton 
(|Biancalana et a/.l . l2QQ4l : ISkrvabin et~al\ , l2QQ3h . 

In a typical supercontinuum setting with /^a > the 
Raman effect increases frequency detuning (A = 5ch—^s) 
between the soliton and its resonance radiation, see 
Fig. EJc). For /^a < the situation is the opposite, 
see Fig. [61(c). The amplitude of the emitted radiation 
is proportional to e~^ 1^1 (a is a constant), where A oc 2 ; 
([Akhmediev and Karlsson', 1995; Biancalana et al., 200; 
[Karpman, 1993; Skryabin et a/., 2003; Wai et al, 19.^ 
Thus for /^a > the soliton self-frequency shift induces 
exponent ial decay of the radiation amplitude with pro 
agation (JBiancalana et all |2004J : iGorbach et al\ , 1200' 
Practically this implies that the soliton emits signifi- 
cant shortwavelength radiation only at the initial stage 
of the supercontinuum generation. The radiation emis- 
sion quickly becomes unnoticeable when the solitons are 
shifted away from the zero GVD point. Then the nat- 
ural question is what causes the continuous blue shift 
of the short-wavelength edge of the supercontinuum, 
see Fig. [3l The answer is - soliton-radiation interac- 



tion ( not mere radiation emission) ([Cumberland et at. , 
I2OO8': 'Gentv et al), '20 04aL [2005[ : 'Gorba ch and Skrvabinl. 
2007b,c; Kudlinski and MussoJ 2008; St one and Knight[ . 
[2008[ : [Travers[ . [2009) . The rest of this section and sections 
V and VI elaborate on this in great details. 



a) 



N 



10 








2nd 
reflection 



± 



1st 
reflection 



200 400 600 

t 




-40 -20 20 40 60 




8 6 4 2 0-2-4-6-8 

Delay (ps) 

FIG. 8 (color online) Experimentally measured XFROG spec- 
trogram showing emission of the strong long wavelength ra- 
diation and its reflection fro m the 2nd soliton in a fiber with 
/^s < (JEfimov et all l2004l ) . The strong resonant radiation 
(wave 1) is emitted by the soliton 1. The weaker soliton 2 re- 
flects a part of this wave backwards and simultaneously trans- 
forms frequency of the reflected radiation (wave 2) . Doted line 
marks the zero GVD wavelength. 



FIG. 7 A weak gaussian pulse in normal GVD range initially 
placed behind the soliton experiences multiple reflections from 
the latter. After each collision some light gets through the 
soliton and some is reflected back. The reflected light is lo- 
calized close to the soliton and almost indistinguishable from 
it in the (z,t)-plot (a). However, the reflected light is clearly 
visible in the spectral (z,^)-plot (b), because its frequency is 
blue-shifted after each collision. 



If Ps > 0, then the resonant radiation emitted by the 
soliton has group velocity less than the soliton itself, so 
that the radiation appears behind the soliton (Fig.[5fa)). 
However, after some propagation the radiation catches 
up with the soliton, which is continuously decelerated by 
the Raman effect, and the two collide (Fig. [5ja), [Tl^a)). 
During this collision the radiation is reflected from the 
soliton backwards, and so the next collision becomes un- 
avoidable, through the same mechanism. The process 
can be repeated many times, see Fig. [71 Note, that the 
radiation colliding with the soliton, should not be nec- 
essarily emitted by the soliton itself, it can be initiated 
by other mechanisms, e.g., via SPM or emitted by other 
solitons (Figs. I7|8p . An important condition for the re- 
flection of the radiation from a soliton to actually happen 
is that the radiation frequency should belong to the nor- 
mal GVD range. The radiation in the anomalous GVD 
range practically does not see solitons, since the model is 
close to the ideal integrable NLS. 

Reflection of the radiation backwards (towards smaller 
z and larger t) from the soliton, implies that the radi- 
ation group velocity is further reduced, i.e. the group 



index dsk for the radiation increases. For the normal 
GVD (9|/c > 0) the increase in the group index happens 
together with the increase in frequency S. Therefore the 
backward reflections have to be accompanied by the blue 
shifting of the radiation frequency, see Fig. [Tfb). This is 
exactly the type of process happening between the soli- 
tons and the short wavelength radi ation in the exp a nding 
supercontinuum shown in Fig . [3] (JGorbach et all l2QQ6l : 
iGorbach and Skrvabinl . l2QQ7bl la). Obviously the reflec- 
tion process is nonlinear in its nature, which is explained 
in details in the next section. 

For /^s < 0, the detuning between the radiation and 
soliton gets smaller with propagation (Fig. [6fc)), there- 
fore the radiation is exponentially amplified in z. This 
amplification leads to the strong spectral recoil on the 
soliton followed by the compensati o n of the soliton self- 
frequ e ncy shift (iBiancalana e^ aU 2004 ;_ Efimov et at. , 



20041: ISkrvabin et all l2003l : iTsov and deSterke CM. . 
2006) (see Fig. [6jb)). In this case the radiation is emit- 
ted ahead of the soliton and has no chance of interacting 
with it, see. Fig. [6l^a). The radiation can, however, in- 
teract with other solitons present in the fibers, bounce 
back from them and has its frequency t ransfo rmed 
(JEfimov et all 120041 : iGorbach and Skrvabinl . l2007c[ ). see 
Fig. El 

We also note at the end of this section that dark soli- 
tons existing for the normal GVD are known to emit 
Chere nkov radiation into the an omalous GVD range 
(Af an as lev et all 1 19961 : iKarpmanl . [l993). Amplification 
of this radiation by the Raman effect and its possible 
use for generation of broad continua has been recently 
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FIG. 9 (color online) Wave number matching diagrams (a,c) 
and XFROG spectrograms (b,d) for the fiber pumped with 
the soliton and cw (Raman effect is disregarded). (a,b) corre- 
spond to /^s > and (c,d) to /^s < 0. The black curved lines 
in (a,c) are k = /32^^/2 + /SsS^ /6. The upper and lower hori- 
zontal lines in (a,c) correspond to the FWM resonances given 
by Eq. (p^ and ([15]), respectively. The middle horizontal 
lines in (a,c) give Cherenkov resonances, Eq. (|lQp . 



investigated (|Milianet fl/.LIl986f ). 



V. SCATTERING OF RADIATION FROM SOLITONS 
AND SUPERCONTINUUM 

In its essence the scattering of a dispersive wave 
from a soliton is a four- wave mixing (FWM) nonlin- 
ear proce ss sensitive to tii e phase matching con d itions 
(lEfimov et all I2QQ4] [2Qq5 I2QQ5I: [C^rbach et all l2QQ6l : 



ISkrvabin and Yulinl . I2QQ5I : lYulin et all l2QQ4l ). The lat- 
ter, however, work out in an unusual way. An impor- 
tant difference of the soliton-radiation interaction with 
the text-book fo ur- wave mixing of dispersive waves only 
(|AgrawallJ2QQ7l ), is that one of the participating fields is 
a non-dispersing pulse (soliton) having straight-line dis- 
persion and moving with a group (not phase) velocity, 
see Eq. (O. 

Lets assume that Fp and dp are the amplitude and 
frequency of the dispersive wave incident on the soliton 
As , while F is the reflected signal fi eld with an unknown 
frequ ency 6. The equation for F is ([Skrvabin and Yulinl , 
l2QQ5h : 
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-l\As\'Fp 

-jA^f; . (11) 



Thus F is exited by \As\ j. p ^^^^ — s-p 



'^Fp and A^F^ with relative effi- 
ciency of the two excitation channels determined by the 
phase matching. Assuming that the incident wave is 
Fp = Cpc'^^^p^-'^p^ and using Eqs. ©-([TD, we find the 
spectral content of the FWM terms: 
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where f{S) = J^^tl;^{t)e'^'dt. 

Using Fourier expansion F = J^ e{6)e^^^^^^~'^^^d8 and 
taking z derivatives of the phases involved (cf. Eq. ([7])) 
we equal the wavenumber of F to the wavenumb ers of 
([T2D a nd (fT3D. The r e sult i s (.Efimov et aL .2001 



2QQ5I: iGorbach et all l2QQ6l : ISkrvabin and Yulinl . 
Yulin et all l2QQ4f ) 



2006 . 



20051 : 



m 
m 



ks{S) 
ks{S) 



[ks{Sp) 
[ks{Sp) 



k{Sp)], 
k{Sp)]. 



(14) 
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Here k is the wavenumber of the generated wave. ks{S) 
and ks{Sp) are the soliton wavenumbers at the generated 
frequency and at the frequency of the wave incident on 
the soliton, respectively. For Ss = 0, ks{S) = q and 
Eqs. (p!4l [T5]) are simplified to the form k{S) = k{Sp), 
k{S) = 2q-k{Sp). 

Solving Eqs. ([T4|) . ([T5]) graphically we find that 
they predict upto 4 res onances, see Fig. [9l^a,c) 
(| Skrvabin and Yulinl . |2QQ5[ ). One resonance is obvious 
6 = 6p and it coincides with the frequency of the incident 
wave (cw pump). |AspFp term and Eq. (p!4j) are respon- 
sible for two nontrivial resonances falling into the regions 
of normal and anomalous GVD. The former one is typ- 
ically much stronger and corresponds to the reflection 
of the wave from the soliton potential. The remaining 
one and the resonance predicted by Eq. (p!5|) usually do 
not scatter much of the incident wave and produce weak, 
but detectable signals. Figs. [9fa,b) show the case when 
the incident wave is reflected backwards from the soliton 
with simultaneous upshift of the reflected wave frequency. 
Figs. [9l(c,d) show the case when the radiation-soliton col- 
lision happens at the front edge of the soliton and the 
wave is reflected ahead of the soliton. The frequency is 
down shifted in the latter example. 

Eqs. ([m [T5|) apply without a change if the soliton 
and dispersive waves are orthogonally polarized, which 
has been used in the experi mental measurements of th e 
soliton-radiation interaction (jEfimov et all l2QQ6l . l2QQ5l ) . 
Thes e measurements fully confirmed validi ty of both Eg . 
(ITil) (JEfimov et all l2QQ5h and Eq. ([15]) (JEfimov et all 
|2QQ6[ ). Corresponding examples of the experimental 
XFROG spectrograms are shown in Fig. ([T 



It has been verified that the frequency up-shift of the 
radiation, resulting from the cascaded back reflection of 
the radiation from intensity of the accelerating soliton 
I Asp, is the mechanism ensuring blue s hift of the short 
wave length edge of the supercontinuum (|Gorbach et all 
|2QQ6[ ). Phase matching conditions (p!4|) work out in a 
way that with every reflection the frequencies of the in- 
cident and reflected waves both tend towards the limit 
point, where the group velocity of the dispersive wave 
coincides with the soliton group velocity (|Gorbach et all 
l2006l : ISkrvabin and Yulinl . 120051) . see Fig. H Indeed, if 
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FIG. 10 (color online) Experimentally measured XFROG 
spectrograms showing radiation-soliton intera ction, (a) cor- 
responds to the resonance g !:iven by Eg . (ITll) (|Efimov et all 
[2005) and (b) to Eq. {[15]) (Efim ov et aa 12006'). (a) shows the 
case of the forward reflected wave, as in Figs. ©(d). Both 
(a) and (b) measurements have been taken for /^s < 0. 



FIG. 11 (color online) Spectral evolution along the fiber 
length. (a)-experiment and (c)-modelling. (b) XFROG 
diagram computed for the propagation distance z = 2m. 
(d) Group index as function of wavelength. The decrease 
in height of the shaded triangles shows the decreasing fre- 
quency difference between the incident and reflected waves 
(^ Gorbach et a/.. . .2006. ). 



the instantaneous soliton frequency is assumed to be zero, 
then the minimum of k{S) corresponds to the frequency 
with the group velocity matched to the soliton group ve- 
locity. Fig. [9fa). If the frequency of the incident wave 
(Sp) is exactly at this minimum or very close to it, then 
the resonance predicted by Eq. (p!4|) is either degenerate 
or close to the degeneracy with 6p. Under such conditions 
the process of frequency conversion practically happens 
within a single wave packet and can be te r med as the 
intrapulse four- wave mixing (|Gorbach et all |2006[ ) . 

In the above picture of the soliton-radiation interac- 
tion, it remains unclear what sustains the efficiency of 
the process over the long propagation distances. Strong 
normal GVD acting on the dispersive wave packet, should 
lead to its sufficiently fast spreading and drop in the peak 
intensity terminating any nonlinear int eraction. In stead, 
as it has been originally measured in (|Hori et all 12004) 
and can be seen from the modelling results (Figs. [3]^b), 
[lUb)), the wavepackets on the short- wavelength edge of 
the continuum remain localised on the femtosecond time 
scale and propagate in the soliton-like regime albeit in 
the normal GVD range. This problem is addressed in 
the next section. 



VI. RADIATION TRAPPING AND SUPERCONTINUUM 

So, what does stop dispersive spreading of the radia- 
tion wavepackets at the short wavelength edge of super- 
continuum and keeps them localised on the femtosecond 
time scale? The short answer is - the refractive index 
change created by the decelerated soliton exerts a special 
type of inertial force and ensures dispersionless propaga- 
tion (trapping) of the radiation. 



Reflection of radiation from a soliton plays the ma- 
jor role in the trapping mechanism. Mathematically 
this is described by Eq. (pTj) with the first nonlin- 
ear term in the right-hand side, see the previous sec- 
tion. Switching into the reference frame moving to- 
gether with the soliton reveals two distinct propagation 
regimes. If the offset of group velocities of the soli- 
ton and radiation is sufficiently large and fiber length 
is relatively short, then the reflection from the soli- 
ton of course happens, but recur rent collisions lead- 
ing to trapping can be disregarded (ISkrvabin and Yulinl , 
|2QQ5[ ). The trapping phenomenon becomes the dom- 
inant feature of the propagation, when group ve- 
locities of the soliton and radiation are sufficiently 
nz — ^ ; ttt; rT~i i ^^^„ n . in rr—— -i r— — — n 



close (Gorba ch and Skrvabin, 2007a. b.c; Hill et a/. , l2QQ9t 



IStone and Kni ght. 2008: Travers and Tavlon |2009[ ). 

Formal consideration of the problem starts from Eq. 
(pT]) complemented by the equation for the soliton field. 
We choose the reference frequencies of the soliton (Ss) 
and radiation (Sp) so that the group velocities are 
matched {k'{Ss) = k'{Sp)) across the zero GVD point 
and assume 

E = As ex.-p[ik{ds)z — idst] + F ex.-p[ik{Sp)z — i6pt]. (16) 

For the amplitudes As(z,t) dind F(z,t) we make substi- 
tutions (|Gorbach and Skrvabinl l2007bl la) 



As = V^(z,Oexp 
F = (/)(2:,0exp 



_,tg + zgz+-^^V 
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_,tp + ,Az + ^^V 



(17) 
(18) 



ts = gz^/2 ,^ = t-ts ,k''= djk{5s^p) . 



We have introduced g parameter to indicate that the 
deceleration rate for the soliton interacting with radia- 
tion can be different from qq for a pure soliton, see Eq. 
(J5j). GVD for the soliton and radiation are anomalous 
ik'^ < 0) and normal {k^ > 0), respectively. Thus direc- 
tions of the frequency shifts acquired by the two waves 
are opposite, see the first terms in the exponential factors 
of Eqs. (p!7|) . (p!8]) . The group indices felt by the soliton 
and radiation {dzts = gz) increase with the same rate 
(see also discussion in the previous section) and hence 
they experience equal negative acceleration. 

Now our problem is reduced to a formal question: 
is there a solution for (j) retaining its localized form? 
In order to answer this formally we substitute ([16]- 
[T8|) into Eq. (1) and neglect several small terms, 
including the ones nonlinear in the ra diation ampli- 
tude l^p (JGorbach and Skrvabinl . l2007c[ ). Then, the z- 
independent, i.e. shape preserving, radiation waves have 
to satisfy 

- \k';dl<Pn + 2j\xP\^4>n + f|^„ = A„0„ ,k';>0. (19) 

The above is the linear Schrodinger equation with the 
effective potential energy V = 27|?/;p + |^. The first 




FIG. 12 (color online) (a) Full line shows the potential V. 
The dashed line shows the potential with the contribution 
from the inert ial force disregarded, which removes the min- 
imum and makes the radiation trapping impossible. Level 
lines indicate ground state, 2nd and 10th modes, (b) shows 
the ground state mode of the potential V c) shows 2nd and 
10th modes. 



FIG. 13 (color online) Spectral evolution of the soliton- 
radiation bound states with n = 1 (a) and n = 3 (b) given by 
Eqs. jni, ([HI). 



term inside V is the repelling potential created by the 
refractive index change induced by the soliton. The dis- 
persive wave reflects from it, as described in the previous 
section. The second term in V is the potential linearly 
increasing in ^, which exists only due to the fact that 
we have switched into the non-inertial frame of reference 
accelerating together with the soliton. Hence this term 
represents a type of inertial force acting on photons. It is 
known from classical mechanics, that inertial forces act 
as usual ones, but show up in the equations of motion 
only, when an appropriate non-inertial frame of reference 
is introduced. 

Overall the potential V(^) has the well defined min- 
imum and therefore supports quasi-localised modes 
(bound states), see Fig. [121 These modes can be 
found either numerically or u sing a variational approach 
([Gorbach and SkrvabinI , |2007c) . Taking soliton plus one 
of these modes and substituting them into Eqs. (J17ll8p 
results in the spectral evolution shown in Fig. [131 
The soliton spectrum moves continuously to the smaller 
frequencies and the radiation spectrum moves towards 
higher frequencies. The sign of /3s is implicit, but very 
important here. If /^a > 0, as in the typical supercntin- 
uum generation experiments, then initially Ss < Sp and 
radiation and soliton spectrally diverge with propagation. 
However, if /3s < 0, then Ss > Sp, while the spectral 
shifts still act in the same directions, so that frequencies 
of the radiation and so l iton co nverge with propagation 
([Gorbach and Skrvabinl . l2007c[ ) . 




(b) 
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FIG. 14 (color online) Spectral (a) and time-domain (b) evo- 
lution of the soliton and trapped radiation modeled using the 
realistic frequency dependence of GVD (cf. (a) with Fig. 
[13]). (b) is plotted in the reference frame accelerating with 
go, which explains the curvature of the soliton trajectory op- 
posite to the one in Fig. [SJa). 



One can notice that taking the higher order modes 
leads to temporal (Fig. [Efc)) and spectral (Fig. [TH]) 
broadening of the radiation. Spectral trajectories in Fig. 
[T3l follow the straight lines because Eqs. (fTT]) . ([TS]) as- 
sume frequency independent GVD. In a real fiber the 
soliton and radiation moving away from the zero GVD 
point encounter increasing absolute values of the GVD. 
This leads to the adiabatic broadening of the soliton and 
reduces its Raman shift (Fig. [HTb)). The result is the 
gradually slowing spectral divergence of the soliton and 
radiation (Fig. [HT a)). Physically, the frequency conver- 
sion of the radiation wavepacket is due to the intrapulse 
four- wave mixing (see previous section), which is made 
possible by the sustained overlap of the radiation and 
soliton pulses. 

Each of the solitons inside the supercontinuum shown 
in Figs. [2fb), [Til has its own radiation pulse contin- 
uously drifting towards shorter wavelengths. We have 
found that the strongest soliton on the long-wavelength 
edge of the supercontinuum spectrum in Fig. 1(a) creates 
the potential V trapping around 20 modes on the short- 
wavelength edge. The radiation captured by the soliton 
can be represented as a superposition of these modes. 
Adiabatic transformation of the soliton power and width 
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with propagation, caused by the increasing dispersion, in- 
duces weak adiabatic evolution of the mode parameters, 
but apart from this the modes are stationary solutions 
and hence their temporal and spectral dynamics are sup- 
pressed. Therefore approximating the radiation field as 



n^,o 
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(20) 



FIG. 15 (color online) (a) Series of images of the condensate 
bouncing off a light sheet (Bongs et aL, 1999). Time interval 
between the images is 2ms along the horizontal axis, (b) 
Numerical modeling of the dispersive radiation bouncing off 
the decelerating soliton. ^ is the dimensinless time measured 
in the reference frame moving together with the soliton and 
z is the distance along the fiber (c) frequency, S, convers ion 
of the bouncing radiation. (|Gorbach and Skrvabinl l2007al l 



gives the very good matching with numerically computed 
spectral evolution of the short wavelength edge of the 
supercontinuum (Gorbach and Skryabin, 2007b). 

At the end of this section we'd like to draw your atten- 
tion to few points. Firstly, the potential barrier, see Fig. 
([12]), on the soliton side is high, but still finite, so that 
some light leaks through it. The leaked radiation is espe- 
cially noticeable in the higher-order modes. Figs. [2fb), 
[T2lf c). Secondly, from numerical modelling of supercon- 
tinuum it can be found that the rate of the soliton self- 
frequency shift on the red edge of the supercontinuum 
is actually higher than for the bare soliton. Analytically 
this effect can be captured if g i s calculated with non- 
linear in ^ terms accounted for ( Gorbach and Skrvabinl , 
l2007cl ). The resulting expression is ^ = go{l-\-Pb'^). Here 
P is the peak power of the radiation and 6^ is a constant. 
Thirdly, to generate spectra shown in Figs. [3land[TlTa). 
(c) we have used the same dispersion profiles and input 
powers, but different input wavelengths. For the Fig. 
[TTJa), (c) the pump was only lOnm away from the zero 
GVD point, this has led to formation of less powerful 
solitons. Hence their frequency shift and associated fre- 
quency shift of the shortwavelength radiation were sub- 
stantially smaller leading to much narrower continua. 
Fourthly, it has been reported that for high pump powers 
the solitons at the infrared edg e of the continuum tend t o 
collide and form bound states (jPodlipenskv et a/.Ll2008[ ). 
These effects are naturally expected to influence radia- 
tion at the short wavelength edge, through the change of 
the trapping potential and the col l ision induced changes 
in soliton frequencies (jLuan et ali uOOw . 



It is important to note, that the simultaneous and op- 
posite soliton-radiation frequency conversion and radi- 
ation trapping have been observed in few experiments 
not related to the mainstream of the fiber supercon- 
tinuum research. First kno wn to us paper is t he 1987 
experiment by Beaud et al (Bea ud et all ll98 7V Then 
there was a gap until 2001, when Nishizawa and Goto 
reported a series of spectral and time domain mea- 
surements of the effect of pulse trappin g by a s oliton 
acros s the zero GVD point (Nishizawa an d Gotj . l200ll 
|2002[ ). For recent experimental observations and fre- 
quency conversion applications of the trapping effect 



see, e.g. jC umberland et all l2QQ8c iHill et al 



Kudlinski and Mussotl. [200^ 
Stone and Knightl . l200g ' 



Nishizawa and Ito 



a 



2009; 



20091 : 



VII. GRAVITY-LIKE EFFECTS, FREAK WAVES AND 
TURBULENCE OF LIGHT IN FIBER 

Parameter g used in the above section has an obvi- 
ous analogy with the acceleration of free fall. Eq. (p!9|) 
can be interpreted as the equation for a quantum par- 
ticle in the gravity field and subject to the additional 
potential created by the soliton. If the radiation wave 
is prepared in such a way that it is both well localised 
in time and shifted away from the potential minimum, 
then it is a highly multimode state. Thus, according to 
the correspondence principle, one should expect quasi- 
classical dynamics to be seen. The wave packet should 
roll down a linear potential towards the soliton, reflect 
back from the latter and, after some time, reconstruct 
itself in the original location. It is known though, that 
the quasi-classical bouncing carries on only for limited 
time, until it is replaced by the complete delocalization 
of the wave packet, which restores its shape a gain later . 
This effect is known as quantum bouncing (|RobinettL 
l2004f ). Previously bouncing has been observed with 
clouds of Bose-condensed ultracold atoms subject to the 
field of grav i ty and reflecting of an atom mirror, see 
(iBongs et all fl999l : ISaba et all Il999[ ) and Fig. Hg^a). 
Fig. [TST b) shows numerically modeled space-time evo- 
lution of the radiation pulse in the normal GVD range 
of an optical fiber bouncing o n a decelerating soliton 
(JGorbach and Skrvabinl l2007al ). Similar light bouncing 
effects have also been reported f or the c urved waveguide 
arrays (JDella Valle et all 120091 : iLonghil . |2008[ ). An im- 
portant feature of our case, is that on each refiection 
from the soliton the frequency of the radiation is up 
shifted (see Fig. [TST c)). Recently, refiection of the radi- 
ation from the soliton and the associated blue shift have 
been interpreted as the frequency shift at the white-hole 
horizon (Philbi n et ali 12008 ). The same work has pre- 
dicted that the quantum effects of horizons, in partic- 
ular Hawking radiation, can potentially be seen due to 
soliton-radiation interaction in optical fibers. 

Another problem recently posed by the supercontin- 
uum and soliton research has been the question about ex - 
istence of optical freak or rogue waves feolli et ali . 120071 ) . 
This phenomenon has been actively studied in the con- 
text of ocean waves, where the rare waves, with proba- 
bility not described by the tails of the gaussian distribu- 
tion, and the amplitude few times larger than the average 
(for the current conditions) wave hight pose serious and 
hardly predictable threat for ships and offshore Indus- 
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tries. NLS model is know n to describe de ep water waves 
including the freak events (|janssenl , |2QQ3[ ) . Therefore one 
can expect appearance of similar phenomena in fiber op- 
tics. In particular, cases of the notable pulse to pulse 
fiuctuations of the supercontinumm, can be attributed to 
generation of the infrared solitons with unusually lar^ e 
amplitudes (JDudlev et all I2QQ8I : ISolli et all l2QQ8l 120071 ) . 
Probability of this to happen is described by the tail 
of the L-shaped distribution function. These freak soli- 
tons emerge essentially due to anomalously strong fo- 
cusing developing in the course of modu lational insta- 
bility and the higher or d er sol iton fission (|Dudlev et all 
l2QQ8l : ISolh et all l2008l . 120071 ). There is also another 
class of localised freak wave solutions of NLS equation. 
These solution are breathers, i.e. localized waves peri- 
odically absorbing and releasi ng their energy in t o the 
continuous wave backg!:round ([Akhmediev et all l2009l : 
lAkhmediev and W abnitz. '1992') . This type of waves has 
the property of sudden appearance and disappearance, 
which is a known feature of the oceanic freak waves, but 
still has not been seen in optical fibers. 

Another area, where analogy between fiuid mechan- 
ics and fiber optics is starting to produce inte r esting 
results is the turbulence resea rch ([Barviau et all l2008l , 
l2009bl : iDvachenko et all Il992[ ) . Spectral broadening due 
to multiple four- wave mixing processes of random weakly 
nonlinear waves can be associated with irreversible evo- 
lution of the spectrally narrow pump towards spec- 
trally broad thermodynamic equilibrium tearviau et all 
[2008, 2009a, b). Current theoretical approaches to the 
turbulence in general and to the turbulent supercon- 
tinuum still have not progressed to the level where 
the role of mixing of incoherent wave fields with co- 
herent wave structures (s olitons) in spectral bro aden- 
ing can be fully revealed ([Pvachenko et all Il992l ). At 
the same time fiber supercontinuum research suggests 
that this difficult case is the most practically rele- 
vant. Complex 'far from equilibrium' dynamics is also 
well known in the fiber based systems, such as fiber 
lasers and coherently pumped fiber resonators, where 
soliton pu lses and spe c tral broadening often coexist , 
see, e.g.. (JBabin et all l2QQ8|; IChouh and Grelul. 12009 : 



1996. : Peng et all 120081 ) 



Kozvre ff e^ a/.L T2OO9I: iLee et all l2005l : iMitschke et~al. . 



VIM. SUMMARY 

For convenience of our readers we summarize here 
those of the soliton properties, which are most impor- 
tant for supercontinumm generation 

• Interaction of a soliton with dispersive radiation 
leads to the phase-matching sensitive generation of 
new frequencies. The most pronounced, out of few 
possible interaction channels, is the refiection of the 
radiation from the refractive index change created 
by the soliton intensity. The refiection happens 
providing the radiation frequency belongs to the 



normal GVD range. Depending on the sign of the 
3rd order dispersion and frequency of the incident 
radiation, the frequency of the refiected wave gets 
either up- or down-shifted. 

• Raman effect decelerates solitons and down-shifts 
their frequency. Such solitons can interact with ra- 
diation repeatedly, trap it on the time scales of 
lOOfs and continuously up-shift the radiation fre- 
quency. 

The prevalent scenario of the supercontinuum generation 
in photonic crystal fibers pumped by femtosecond pulses 
with the input wavelength around the zero GVD point 
can be summarized as: 

• Spectrum of the input pulse is distributed in some 
proportion between the frequency ranges with nor- 
mal and anomalous GVD. This happens through 
the combination of nonlinear processes. SPM dom- 
inates during the first few centimeters of propaga- 
tion. Then the soliton fission accompanied by the 
radiation emission and refiection of the dispersive 
waves from emerging solitons lead to further spec- 
tral broadening. 

• The next stage is when the Raman shifted solitons 
on the long- wavelength edge of the supercontinuum 
enter into the regime of the cascaded interaction 
with dispersive radiation. This quickly leads to 
formation of the bound soliton-radiation states re- 
sponsible for continuous spectral divergence of the 
supercontinuum edges. The necessary condition for 
this to happen is the near matching of the group 
velocities across the zero GVD point. 



Using nano-second or cw pump for supercontinuum 
generation leads to modulational instability and sub- 
sequent creation of a soliton train. The latter traps 
radiation and the above scenario is realized again al- 
beit with greater number of s o litons ([Cumberland et al . 
20081: iKudlinski and MussotL 120081 : IStone and KnightL 
20081: iT ravers'. '200 9|) and greater sensitivi ty to noise 
(jSolH et al . 2007 ; Turitsvn and Derevvanko. |2008[ ). 



IX. PERSPECTIVES 

Interaction between solitons and radiation, optical tur- 
bulence, freak waves, and development of ideas around 
the gravity-like forces exerted on light by solitons all 
are on the Hst of problems stimulated by the fiber su- 
percontinuum research and undergoing the stage of ac- 
tive exploration. Supercontinuum generation has been 
of course known outside the fibe r context in bulk 
solids, liquids and gases, s ee, e. g . (iBerge et all 12007 ; 
Couairon and Mvsvrowiczl 120071 : iKolesik and Molonev , 
2008[ ). where all three space dimensions are impor- 
tant. Role of spatial and spatio-temporal solitons 
in these systems is still far from been fully explored 
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(JBerge et all I2QQ7I : lYulin et all |2QQ5[ ). Generation of 
broad spatial and frequency spectra in nonlinear pho- 
tonic crystals is another area where interaction of 
solitons with diffracting and dispersing waves c an be 
important (Babushkin et al, 2007: iBartal et al , 2006 : 



Benton et all 12008^ iDong et all 120081 : IJia et all boO? : 
Manela et all |2006[ ). Strong field localisation in metalic 
optical nano-antennas has be en demonstrated to lead t o 
supercontinuum generation ([Muhlschlegel et all . |2005[ ), 
thereby linking the effects discussed above with nanopho- 
tonics. A possibility of nonlinear optical processes in 
nature made optical waveguides, e.g., found in sea or- 
ganisms and sometimes having a structure of fibers with 
few micron core diameters and photonic cry stal cladding 
(Aizenberg e t all I2004J : iKulchin et all 1200 8") , remains an 
intriguing problem to consider. Overall, a bidirectional 
fiow of ideas between fiber photonics and other branches 
of optics and physics in general is more than likely to 
stimulate further progress in the soliton and supercon- 
tinuum related research. 



References 

Afanasjev V.V., Y.S. Kivshar, and C.R. Menyuk, Opt. Lett. 

21, 1975 (1996). 
Agrawal, G., 2007, Nonlinear Fiber Optics (Academic Press). 
Aizenberg, J., V. C. Sundar, A. D. Yablon, J. C. Weaver, and 

G. Chen, 2004, Proc. Nat. Acad. Sciences of the USA 101, 

3358. 
Akhmediev, N., A. Ankiewicz, and M. Taki, 2009, Phys. Lett. 

A 373, 675. 
Akhmediev, N., and M. Karlsson, 1995, Phys. Rev. A 51, 

2602. 
Akhmediev, N., W. Krolikowski, and A. Lowery, 1996, Opt. 

Commun. 131, 260. 
Akhmediev, N. N., and S. Wabnitz, 1992, J. Opt. Soc. Am. 

B 9, 236. 
Amiranashvili, S., A. G. Vladimirov, and U. Bandelow, 2008, 

Phys. Rev. A 77, 063821. 
Babin, S. A., V. Karalekas, E. V. Podivilov, V. K. Mezentsev, 

P. Harper, J. D. Ania- Cast anon, and S. K. Turitsyn, 2008, 

Phys. Rev. A 77, 033803. 
Babushkin, L, A. Husakou, J. Herrmann, and Y. S. Kivshar, 

2007, Opt. Express 15, 11978. 
Bartal, G., O. Manela, and M. Segev, 2006, Phys. Rev. Lett. 

97, 073906. 
Barviau, B., B. Kibler, S. Coen, and A. Picozzi, 2008, Opt. 

Lett. 33, 2833. 
Barviau, B., B. Kibler, A. Kudlinski, A. Mussot, G. Millot, 

and A. Picozzi, 2009a, Opt. Express 17, 7392. 
Barviau, B., B. Kibler, and A. Picozzi, 2009b, Phys. Rev. A 

79, 063840. 
Beaud, P., W. Hodel, B. Zysset, and H. Weber, 1987, IEEE 

J. Quant. Elect. 23, 1938. 
Benton, C. J., A. V. Gorbach, and D. V. Skryabin, 2008, 

Phys. Rev. A 78, 033818. 
Berge, L., S. Skupin, R. Nuter, J. Kasparian, and J. P. Wolf, 

2007, Rep. Prog. Phys. 70, 1633. 
Biancalana, P., D. V. Skryabin, and A. V. Yulin, 2004, Phys. 

Rev. E 70, 016615. 



Bongs, K., S. Burger, G. Birkl, K. Sengstock, W. Ertmer, 

K. Rzazewski, A. Sanpera, and M. Lewenstein, 1999, Phys. 

Rev. Lett. 83, 3577. 
Chouli, S., and P. Grelu, 2009, Opt. Express 17, 11776. 
Corwin, K. L., N. R. Newbury J. M. Dudley S. Coen, S. A. 

Diddams, K. Weber, and R. S. Windeler, 2003, Phys. Rev. 

Lett. 90, 113904. 
Couairon, A., and A. Mysyrowicz, 2007, Phys. Rep. 441, 47. 
Cristiani, L, R. Tediosi, L. Tartara, and V. Degiorgio, 2004, 

Opt. Express 12, 124. 
Cumberland, B. A., J. C. Travers, S. V. Popov, and J. R. 

Taylor, 2008, Opt. Lett. 33, 2122. 
Delia Valle, G., M. Savoini, M. Ornigotti, P. Laporta, 

V. Foglietti, M. Finazzi, L. Duo, and S. Longhi, 2009, Phys. 

Rev. Lett. 102, 180402. 
Dong, R., C. E. Ruter, D. Kip, O. Manela, M. Segev, C. L. 

Yang, and J. J. Xu, 2008, Phys. Rev. Lett. 101, 183903. 
Dudley J. M., and S. Coen, 2002, Opt. Lett. 27, 1180. 
Dudley J. M., G. Genty and S. Coen, 2006, Rev. Mod. Phys. 

78, 1135. 
Dudley, J. M., G. Genty, and B. J. Eggleton, 2008, Opt. Ex- 
press 16, 3644. 
Dudley, J. M., and J. R. Taylor, 2009, Nature Photonics 3, 

85. 
Dyachenko, S., A. C. Newell, A. Pushkarev, and V. E. Za- 

kharov, 1992, Physica D 57, 96. 
Efimov, A., A. J. Taylor, F. G. Omenetto, A. V. Yulin, N. Y. 

Joly, F. Biancalana, D. V. Skryabin, J. C. Knight, and P. S. 

Russell, 2004, Opt. Express 12, 6498. 
Efimov, A., A. J. Taylor, A. V. Yulin, D. V. Skryabin, and 

J. C. Knight, 2006, Opt. Lett. 31, 1624. 
Efimov, A., A. V. Yulin, D. V. Skryabin, J. C. Knight, N. Joly 

F. G. Omenetto, A. J. Taylor, and P. Russell, 2005, Phys. 

Rev. Lett. 95, 213902. 
Falk, P., M. H. Frosz, and O. Bang, 2005, Opt. Express 13, 

7535. 
Frosz, M. H., P. Falk, and O. Bang, 2005, Opt. Express 13, 

6181. 
Gagnon, L., and P. A. Belanger, 1990, Opt. Lett. 15, 466. 
Genty, G., M. Lehtonen, and H. Ludvigsen, 2004a, Optics 

Express 12, 4614. 
Genty, G., M. Lehtonen, and H. Ludvigsen, 2005, Opt. Lett. 

30, 756. 
Genty, G., M. Lehtonen, H. Ludvigsen, and M. Kaivola, 

2004b, Opt. Express 12, 3471. 
Gorbach, A., D. Skryabin, J. Stone, and J. Knight, 2006, Opt. 

Express 14, 9854. 
Gorbach, A. V., and D. V. Skryabin, 2007a, Opt. Express 15, 

14560. 
Gorbach, A. V., and D. V. Skryabin, 2007b, Nature Photonics 

1, 653. 
Gorbach, A. V., and D. V. Skryabin, 2007c, Phys. Rev. A 76, 

053803. 
Gorbach, A. V., and D. V. Skryabin, 2008, Opt. Express 16, 

4858. 
Gorshkov, K. A., and L. A. Ostrovsky, 1981, Physica D 3, 

428. 
Gu, X., M. Kimmel, A. P. Shreenath, R. Trebino, J. M. Dud- 
ley S. Coen, and R. S. Windeler, 2003, Opt. Express 11, 

2697. 
Harbold, J. M., F. O. Ilday F. W. Wise, T. A. Birks, W. J. 

Wadsworth, and Z. Chen, 2002, Opt. Lett. 27, 1558. 
Hasegawa, A., and F. Tappert, 1973, Appl. Phys. Lett. 23, 

142. 



13 



Herrmann, J., U. Griebner, N. Zhavoronkov, A. Husakou, 

D. Nickel, J. C. Knight, W. J. Wadsworth, P. S. J. Russell, 

and G. Korn, 2002, Phys. Rev. Lett. 88, 173901. 
Hill, S., C. E. Kuklewicz, U. Leonhardt, and F. Konig, 2009, 

Opt. Express 17, 13588. 
Hori, T., N. Nishizawa, T. Goto, and M. Yoshida, 2004, J. 

Opt. Soc. Am. B 21, 1969. 
Husakou, A. V., and J. Herrmann, 2001, Phys. Rev. Lett. 87, 

203901. 
Janssen, P. A. E. M., 2003, Journal Of Physical Oceanography 

33, 863. 
Jia, S., W. Wan, and J. W. Fleischer, 2007, Phys. Rev. Lett. 

99, 223901. 
Karpman, V. L, 1993, Phys. Rev. E 47, 2073. 
Karpman, V.L, Phys. Lett. A 181, 211 (1993). 
Kaup, D. J., and A. C. Newell, 1978, Proc. Royal Soc. Of 

London A 361, 413. 
Kivshar, Y.S., and B.A. Malomed, 1989, Rev. Mod. Phys. 61, 

763. 
Knight, J.C., and D.V. Skryabin, Opt. Express 15, 15365 

(2007). 
Kolesik, M., and J. V. Moloney, 2008, Opt. Express 16, 2971. 
Kozyreff, G., M. Tlidi, A. Mussot, E. Louvergneaux, M. Taki, 

and A. G. Vladimirov, 2009, Phys. Rev. Lett. 102, 043905. 
Kudlinski, A., and A. Mussot, 2008, Opt. Lett. 33, 2407. 
Kulchin, Y. N., S. N. Bagaev, O. A. Bukin, S. S. Voznesen- 

skii, A. L. Drozdov, Y. A. Zinin, L G. Nagornyi, E. V. 

Pestryakov, and V. L Trunov, 2008, Techn. Phys. Lett. 34, 

633. 
Kuznetsov, E. A., A. V. Mikhailov, and L A. Shimokhin, 1995, 

Physica D 87, 201. 
Lee, J., Y. Takushinia, and K. Kikuchi, 2005, Opt. Lett. 30, 

2599. 
Longhi, S., 2008, Phys. Rev. A 77, 035802. 
Luan, F., D. Skryabin, A. Yulin, and J. Knight, 2006, Opt. 

Express 14, 9844. 
Manela, O., G. Bartal, M. Segev, and H. Buljan, 2006, Opt. 

Lett. 31, 2320. 
Milian, C., D.V. Skryabin, and A. Ferrando, Opt. Lett. 34, 

2096 (2009). 
Mitschke, F., and L. Mollenauer, 1986, Opt. Lett. 11, 659. 
Mitschke, F., G. Steinmeyer, and A. Schwache, 1996, Physica 

D 96, 251. 
Mollenauer, L., and J. Gordon, 2006, Solitons in Optical 

Fibers: Fundamentals and Applications (Academic Press). 
Mollenauer, L. F., R. H. Stolen, and J. P. Gordon, 1980, Phys. 

Rev. Lett. 45, 1095. 
Muhlschlegel, P., H. J. Eisler, O. J. F. Martin, B. Hecht, and 



D. W. Pohl, 2005, Science 308, 1607. 
Nishizawa, N., and T. Goto, 2001, Opt. Express 8, 328. 
Nishizawa, N., and T. Goto, 2002, Opt. Express 10, 1151. 
Nishizawa, N., and K. Itoh, 2009, Appl. Phys. Express 2, 

062501. 
Peng, J., F. Zhu, and A. V. Sokolov, 2008, Opt. Lett. 33, 

1620. 
Philbin, T. G., C. Kuklewicz, S. Robertson, S. Hill, F. Konig, 

and U. Leonhardt, 2008, Science 319, 1367. 
Podlipensky, A., P. Szarniak, N. Y. Joly, and P. S. J. Russell, 

2008, J. Opt. Soc. Am. B 25, 2049. 
Ranka, J. K., R. S. Windeler, and A. J. Stentz, 2000, Opt. 

Lett. 25, 25. 
Robinett, R. W., 2004, Phys. Rep. 392, 1. 
Russell, P. S. J., 2006, J. Lightwave Tech. 24, 4729. 
Saba, C. V., P. A. Barton, M. G. Boshier, L G. Hughes, 

P. Rosenbusch, B. E. Sauer, and E. A. Hinds, 1999, Phys. 

Rev. Lett. 82, 468. 
Scott, A., 1999, Nonlinear Science: Emergence and Dynamics 

of Coherent Structures (Oxford). 
Skryabin, D. V., F. Luan, J. C. Knight, and P. S. Russell, 

2003, Science 301, 1705. 
Skryabin, D. V., and A. V. Yulin, 2005, Phys. Rev. E 72, 

016619. 
Solh, D. R., C. Ropers, and B. Jalah, 2008, Phys. Rev. Lett. 

101, 233902. 
Solh, D. R., C. Ropers, P. Koonath, and B. Jalah, 2007, Na- 
ture 450, 1054. 
Stone, J. M., and J. C. Knight, 2008, Opt. Express 16, 2670. 
Travers, J., 2009, Optics Express 17, 1502. 
Travers, J. C., and J. R. Taylor, 2009, Opt. Lett. 34, 115. 
Tsoy, E., and deSterke CM., 2006, J. Opt. Soc. Am. B 23, 

2425. 
Turitsyn, S. K., and S. A. Derevyanko, 2008, Phys. Rev. A 

78, 063819. 
Wadsworth, W. J., A. Ortigosa-Blanch, J. C. Knight, T. A. 

Birks, T. P. M. Man, and P. S. Russeh, 2002, J. Opt. Soc. 

Am. B 19, 2148. 
Wai, P. K. A., H. H. Chen, and Y. C. Lee, 1990, Phys. Rev. 

A 41, 426. 
Yulin, A. v., D. V. Skryabin, and P. S. J. Russeh, 2004, Opt. 

Lett. 29, 2411. 
Yulin, A. v., D. V. Skryabin, and P. S. J. Russeh, 2005, Opt. 

Lett. 30, 525. 
Zakharov, V. E., and A. B. Shabat, 1972, Sov. Phys. JETP 

34, 62. 



